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Abstract 


This  study  was  undertaken  to  numerically  investigate  the  behaviour  of  the  submarine  nickel 
aluminum  bronze  (NAB)  valves  (SW403  and  SW226)  with  corrosion  damage  under  hydrostatic 
loading  at  Defence  R&D  Canada  -  Atlantic  (DRDC  Atlantic),  and  to  validate  the  numerical 
models  using  experimental  results.  Finite  element  analysis,  involving  nonlinear  elastic  plastic 
static  analysis  was  performed  using  the  implicit  solver  in  LS-DYNA  software  system.  The 
experimental  testing  showed  that  the  smaller  SW403  valve  generally  responded  with  lower  strains 
than  the  larger  SW226  valve  at  the  same  load.  In  addition,  the  strains  in  the  body  increased  with 
corrosion  damage,  and  were  much  higher  than  those  in  the  neck.  This  was  expected  since  there 
was  no  corrosion  damage  in  the  neck  region,  there  should  be  no  effect  of  corrosion  damage  on  the 
circumferential  strain  in  the  neck.  The  finite  element  results  for  the  axial  strain  in  the  body  in 
general  showed  reasonably  good  agreement  with  the  experimental  results.  However,  the 
circumferential  strain  in  both  the  body  and  the  neck  were  in  general  underestimated.  This 
indicates  that  further  refinement  of  the  finite  element  mesh  is  required. 


Resume 


La  presente  etude  a  etc  menee  par  R&D  pour  la  Defense  Canada  -  Atlantique  (RDDC 
Atlantique)  afin  d’examiner,  a  I’aide  d’une  simulation  numerique,  la  reaction  des  vannes  en 
bronze  au  nickel-aluminium  (BNA)  (SW403  et  SW226)  endommagees  par  la  corrosion  lors  d’un 
chargement  hydrostatique  et  pour  valider  les  modeles  numeriques  a  I’aide  de  resultats 
experimentaux.  La  methode  d’analyse  par  elements  finis  qui  comprend  une  analyse  statique  non- 
lineaire  elastoplastique  a  etc  effectuee  a  I’aide  du  resolveur  implicite  du  logiciel  LS-DYNA.  Les 
essais  ont  demontre  que,  de  fagon  generale,  la  vanne  de  plus  petite  taille  SW403  presente  moins 
de  deformations  pour  le  meme  chargement  que  la  vanne  de  plus  grande  taille  SW226.  Par  ailleurs, 
le  nombre  de  deformations  dans  le  corps  augmente  avec  I’etendue  de  la  corrosion  et  dies  y  sont 
beaucoup  plus  nombreuses  que  dans  le  col,  chose  qui  etait  a  prevoir.  En  effet,  s’il  n’y  a  pas  de 
dommages  dus  a  la  corrosion  dans  la  region  du  col,  il  ne  devrait  pas  avoir  de  deformations 
longitudinales  dus  a  la  corrosion  sur  le  col.  En  generale,  les  resultats  des  analyses  par  elements 
finis  pour  la  deformation  longitudinale  sur  le  corps  sont  semblablement  les  memes  que  lors  des 
essais.  Cependant,  les  deformations  transversales  sur  le  corps  et  le  col  de  la  vanne  ont  etc  sous- 
estimes,  indiquant  qu’il  serait  necessaire  d’ameliorer  le  reseau  maille  des  elements  finis. 
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Executive  summary 


Hydrostatic  Response  of  Submarine  Nickel  Aluminum  Bronze 
Valves  with  Corrosion  Damage 

B.K.C.  Yuen;  T.S.  Koko;  R.  Warner;  DRDC  Atlantic  CR  2008-105;  Defence  R&D 
Canada  -  Atlantic.  July  2008. 

Introduction  or  background:  This  study  was  undertaken  in  support  of  on-going  mid-life 
evaluations  of  the  Victoria  Class  submarine  nickel  aluminum  bronze  valves.  In  order  for  the 
valves  to  perform  their  intended  functions,  it  is  important  that  the  integrity  (structural  and  water 
tightness)  of  the  valves  is  maintained  under  all  operational  loading  conditions  throughout  the  life 
of  the  vessel.  The  objective  of  this  study  is  to  numerically  investigate  the  behaviour  of  the 
submarine  nickel  aluminum  bronze  (NAB)  valves  (SW403  and  SW226)  with  corrosion  damage 
under  hydrostatic  loading,  and  to  validate  the  numerical  models  using  experimental  results. 
Corrosion  damage  was  simulated  by  gradually  removing  material  from  the  inside  the  valve  body. 
Finite  element  models  of  the  valve  systems  were  then  developed  to  simulate  the  hydrostatic  test 
configuration  developed  at  DRDC  Atlantic.  Finite  element  analysis,  involving  nonlinear  elastic 
plastic  static  analysis  was  performed  using  the  implicit  solver  in  LS-DYNA  software  system. 

Results:  The  experimental  testing  showed  that  the  smaller  SW403  valve  generally  responded 
with  lower  strains  than  the  larger  SW226  valve  at  the  same  load.  In  addition,  the  strains  in  the 
body  increased  with  corrosion  damage,  and  were  much  higher  than  those  in  the  neck.  This  was 
expected  since  there  was  no  corrosion  damage  in  the  neck  region,  there  should  be  no  effect  of 
corrosion  damage  on  the  circumferential  strain  in  the  neck.  The  finite  element  results  for  the 
axial  strain  in  the  body  in  general  showed  reasonably  good  agreement  with  the  experimental 
results.  However,  the  circumferential  strain  in  both  the  body  and  the  neck  were  in  general 
underestimated. 

Significance:  The  numerical  modeling  capabilities  of  the  two  valve  systems  have  been  developed 
to  predict  the  response  of  the  damage  containing  valves  under  hydrostatic  loading  conditions.  The 
numerical  modeling  results  have  shown  reasonably  good  results  in  predicting  the  axial  strains 
under  various  loading  conditions.  Once  properly  validated,  this  numerical  modeling  approach  can 
be  used  for  fitness-for-service  assessments  of  damage  containing  (NAB)  valves  in  the  Victoria 
class  submarine  seawater  systems,  assessment  of  the  current  standards  that  govern  valve 
management  practice,  to  quantify  the  levels  of  conservatism,  and  for  optimization  of  experimental 
valve  assessment  trials. 

Future  plans:  The  effects  of  issues  such  as  mesh  size,  element  formulation,  element  order  and 
material  properties  will  be  investigated.  Addressing  these  issues  would  improve  the  accuracy  of 
the  FE  models  so  that  they  more  closely  match  experimental  measurements  and  further  increase 
confidence  in  their  reliability.  Hydrostatic  testing  of  more  valve  systems  (at  least  two  more)  will 
be  performed.  Numerical  simulations  of  these  valves  will  also  be  performed  to  further  validate 
the  numerical  approach.  Once  validated,  the  numerical  models  can  be  used  to  generate 
performance  charts  for  existing  and  new  valve  designs.  In  addition,  numerical  modelling 
capabilities  for  prediction  of  dynamic  loading  response  of  the  damage  containing  valves  will  also 
be  developed  and  validated  using  dynamic  loading  test  results. 
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Reaction  hydrostatique  des  vannes  de  sous-marins  en  bronze 
au  nickel-aluminium  endommagees  par  la  corrosion 

B.K.C.  Yuen;  T.S.  Koko;  R.  Warner;  DRDC  Atlantic  CR  2008-105;  R  &  D  pour  la 
Defense  Canada  -  Atlantique.  Juillet  2008. 

Introduction  on  contexte  :  La  presente  etude  a  ete  menee  en  soutien  aux  evaluations  continues 
de  mi-duree  des  vannes  en  bronze  au  nickel-aluminium  (BNA)  des  sous-marins  de  classe 
Victoria.  Pour  que  les  vannes  puissent  effectuer  leurs  functions  prevues,  il  est  important  que  leur 
integrite  (structurale  et  leur  etancheite  a  I’eau)  soit  maintenue  dans  toutes  les  conditions  de 
chargement  operationnelles  pour  la  duree  complete  du  sous-marin.  L’objectif  de  la  presente  etude 
est  d’examiner,  a  I’aide  d’une  simulation  numerique,  la  reaction  des  vannes  en  bronze  au  nickel- 
aluminium  (SW403  et  SW226)  endommagees  par  la  corrosion  lors  d’un  chargement 
hydrostatique  et  pour  valider  les  modeles  numeriques  a  I’aide  de  resultats  experimentaux.  Les 
dommages  dus  a  la  corrosion  ont  etc  simules  en  enlevant  graduellement  du  materiau  a  I’interieur 
du  corps  de  la  vanne.  Des  modeles  a  elements  finis  ont  ensuite  etc  elabores  pour  simuler  I’essai 
hydrostatique  congu  par  RDDC  Atlantique.  La  methode  d’ analyse  par  elements  finis,  qui 
comprend  une  analyse  statique  non-lineaire  elastoplastique,  a  etc  effectuee  a  I’aide  du  resolveur 
implicite  du  logiciel  LS-DYNA. 

Resultats  :  Les  essais  ont  demontre  que,  de  fagon  generate,  la  vanne  de  plus  petite  taille  SW403 
presente  moins  de  deformations  pour  le  meme  chargement  que  la  vanne  de  plus  grande  taille 
SW226.  Par  ailleurs,  le  nombre  de  deformations  sur  le  corps  augmente  avec  I’etendue  de  la 
corrosion  et  dies  y  sont  beaucoup  plus  nombreuses  que  sur  le  col,  chose  qui  etait  a  prevoir.  En 
effet,  s’il  n’y  a  pas  de  dommages  dus  a  la  corrosion  dans  la  region  du  col,  il  ne  devrait  pas  avoir 
de  deformations  longitudinales  dues  a  la  corrosion  sur  le  col.  En  generate,  les  resultats  des 
analyses  par  elements  finis  pour  la  deformation  longitudinale  sur  le  corps  sont  semblablement  les 
memes  que  lors  des  essais.  Cependant,  les  deformations  transversales  sur  le  corps  et  le  col  de  la 
vanne  ont  etc  sous-estimees. 

Importance  :  Les  fonctions  de  modelisation  numerique  des  deux  vannes  ont  etc  elaborees  pour 
predire  la  reaction  des  vannes  endommagees  par  la  corrosion  a  un  chargement  hydrostatique.  La 
modelisation  numerique  s’ est  averee  plutot  bonne  pour  predire  les  deformations  longitudinales 
sous  divers  chargements.  Une  fois  bien  validee,  la  modelisation  numerique  pourra  etre  utilisee 
pour  evaluer  I’aptitude  au  service  des  vannes  BNA  endommagees  dans  les  systemes  d’eau  de  mer 
des  sous-marins  de  classe  Victoria,  pour  evaluer  les  normes  actuelles  qui  gouvernent  la  gestion 
des  vannes,  pour  definir  de  fagon  quantitative  les  degres  de  moderation  a  appliquer  et  pour 
optimiser  les  essais  experimentaux  des  vannes. 

Perspectives  :  Les  impacts  des  facteurs  comme  la  dimension  du  reseau  maille,  la  formulation  des 
elements,  I’ordre  des  elements  et  les  proprietes  des  materiaux  seront  evaluees.  L’examen  de  ces 
facteurs  permettrait  d’ameliorer  I’exactitude  des  modeles  a  elements  finis  pour  qu’ils  puissent 
correspondre  davantage  aux  mesures  experimentales  et  permettrait  aussi  d’augmenter  la 
confiance  en  leur  fiabilite.  Au  moins  deux  autres  vannes  seront  soumises  a  des  essais 
hydrostatiques.  Des  simulations  numeriques  des  vannes  seront  egalement  effectuees  pour  valider 


IV 


DRDC  Atlantic  CR  2008-105 


I’approche  numerique.  Une  fois  valides,  les  modeles  numeriques  peuvent  etre  utilises  pour  creer 
des  graphiques  de  performance  pour  les  vannes  actuelles  et  futures.  De  plus,  les  fonctions  de  la 
modelisation  numerique  permettant  de  predire  la  reaction  des  vannes  endommagees  au 
chargement  dynamique  seront  egalement  elaborees  et  validees  a  I’aide  des  resultats  de  I’essai  du 
chargement  dynamique. 
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1  Introduction 


1.1  Background 

As  part  of  the  efforts  to  assess  the  structural  integrity  on  the  submarine  NAB  valves,  DRDC 
Atlantic  has  performed  experimental  and  numerical  studies  to  assess  the  response  of  the  valves 
under  hydrostatic  loading.  In  the  previous  numerical  simulations  [1],  a  number  of  simplifying 
assumptions  and  simplifications  were  made  to  obtain  approximations  of  the  hydrostatic  responses 
of  two  valves,  namely  SW403  and  SW226.  The  previous  experimental  investigations  tested  the 
valves  to  failure  (leakage)  under  hydrostatic  loading,  but  did  not  provide  load-strain  data,  which 
made  it  difficult  to  use  the  experimental  results  to  validate  the  numerical  simulations.  DRDC 
Atlantic  has  now  performed  additional  tests  to  obtain  load-strain  data  under  hydrostatic  loading 
for  a  number  of  NAB  valves.  It  is  therefore  necessary  to  refine  the  numerical  model  and  validate 
them  against  the  experimental  results.  Furthermore,  the  previous  numerical  model  did  not 
account  for  the  presence  of  corrosion  damage  and  it  is  required  to  incorporate  these  into  the 
model  to  assist  in  making  fitness-for-service  decisions.  This  report  describes  the  numerical 
investigations  performed  in  support  of  the  experiments.  The  specific  tasks  to  achieve  these 
objectives  are  listed  below. 

1.2  Objectives  and  Scope 

The  main  objectives  of  this  study  were  to  numerically  investigate  the  behaviour  of  the  submarine 
NAB  valves  with  corrosion  damage  under  hydrostatic  loading  conditions,  and  to  validate  the 
numerical  models  using  experimental  results.  The  scope  of  the  study  includes  the  following: 

1 .  Model  refinement 

This  task  involved  the  refinement  of  the  previous  numerical  models  of  the  SW403  and 
SW226  valves,  in  order  to  remove  some  of  the  simplifications  made  in  the  previous  study  [1]. 
In  that  study,  the  bonnet  was  assumed  to  be  connected  to  the  valve  body  at  the  interface.  This 
reasonable  assumption  was  made  in  order  to  simplify  the  model,  but  it  has  the  effect  of 
making  the  bonnet-valve  interface  somewhat  stiffer  than  the  physical  problem.  A  more 
accurate  representation  would  be  to  model  the  holes  and  bolts  in  detail.  In  the  present  study, 
the  model  was  refined  to  include  detailed  geometry  of  the  bonnet  and  the  bolt  assembly  to 
accurately  model  the  strain  in  the  bonnet.  Using  the  refined  finite  element  models,  analyses 
of  the  hydrostatic  responses  of  the  SW403  and  SW226  valves  were  preformed  and  compared 
to  the  results  of  the  previous  study  and  experimental  results.  Based  on  the  results,  the 
numerical  modeling  strategy  for  hydrostatic  response  analysis  was  prescribed. 

2.  Valve  failure  analysis 

The  following  activities  were  performed  in  this  task: 

a.  Thickness  measurements  of  two  valves  (SW403  and  SW226)  after  some  material 
removal  from  the  inside  the  valve  body  (one  side  only)  were  performed  using  the 
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CMM  machine  at  Dalhousie  University.  The  previous  valve  geometry  data  for  each 
valve  were  used  as  baseline  data  for  the  measurement. 

b.  Development  of  numerical  models  of  two  valves  (SW403  and  SW226)  based  on  the 
new  thickness  data  in  the  predefined  thickness  reduction  region. 

c.  Analysis  of  the  hydrostatic  response  of  the  valves  with  reduced  wall  thickness  at  the 
specified  region.  Comparison  with  the  strain  data  obtained  from  the  experiments. 

1.3  Organization  of  this  Document 

The  remainder  of  this  document  is  organized  as  follows: 

•  Chapter  2  provides  a  description  of  the  problem,  where  details  of  the  NAB  valve 
configurations,  loading  conditions  and  analysis  requirements  are  presented; 

•  Chapter  3  describes  of  the  measurements  taken  to  completely  define  the  valve  geometries 
and  CAD  drawings  of  the  valves  are  presented; 

•  Chapter  4  provides  details  of  the  finite  element  models  used  to  analyze  the  NAB  valves; 

•  Chapter  5  discusses  the  analysis  of  the  NAB  valves  subject  to  hydrostatic  loading;  and 

•  Chapter  6  provides  a  summary  of  the  work,  conclusions  reached  and  recommendations  for 
future  work. 
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2  Problem  Description 


2.1  Configuration  of  Submarine  NAB  Valves 

Figure  1  shows  pictures  of  typical  valves  on  the  Victoria  Class  submarines  and  Figure  2  shows 
pictures  of  the  valves  after  some  material  removal  from  the  inside  the  valve  body  in  this  study. 
These  valves  are  cast  from  nickel  aluminum  bronze  (NAB)  material,  and  are  designed  for  a  safe 
commission  life  (SCL)  of  15  years  [2].  As  the  valves  approach  their  design  SCL,  DRDC  has 
initiated  studies  to  conduct  the  fitness-for-purpose  evaluations  under  various  loading  and 
environmental  conditions,  and  maintenance  and  repair  scenarios.  One  of  these  studies  involved 
experimental  investigations  into  the  behaviour  of  valves  with  corrosion  damage  under  hydrostatic 
loading  conditions.  The  present  study  seeks  to  support  the  experimental  studies,  using  numerical 
simulation.  Description  of  the  experimental  set-ups  is  provided  below. 

2.2  Hydrostatic  Load  Tests 

Figure  3  shows  the  experimental  setup  during  hydrostatic  testing.  Each  valve  was  pressurized 
gradually  and  manually  inside  the  test  container  to  a  preset  pressure  and  then  gradually  decreased 
to  zero.  Each  valve  was  also  instrumented  with  three  strain  gauges  to  capture  the  axial  and  the 
circumferential  strain  in  the  valve  body  and  the  circumferential  strain  in  the  neck,  as  shown  in 
Figure  4.  The  data  on  both  the  hydrostatic  pressure  and  the  strains  were  recorded  during  the  test 
using  a  data  acquisition  system.  The  valves  tested  included  valves  designated  as  SW403  and 
SW226.  In  addition,  the  wall  thickness  of  the  thinner  portion  of  the  valves  was  gradually  reduced 
by  material  removal  as  shown  in  Figure  2,  resulting  in  a  total  of  two  and  four  cases  being  tested 
for  the  SW403  and  the  SW226  valve,  respectively. 

2.3  Requirements 

It  was  required  to  numerically  simulate  the  hydrostatic  responses  of  two  valves  designated  as 
SW403  and  SW226  with  reduced  wall  thickness  at  the  specified  region,  and  to  compare  the 
numerical  results  to  experimental  results,  where  available. 
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Figure  1:  Typical  NAB  Valves  on  Victoria  Class  Submarine. 
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(a) 


(b) 


Figure  2:  NAB  Valves  after  Some  Material  Removal  from  the  Inside  the  Valve  Body:  (a) 

SW403,  (b)SW226  [3]. 
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Figure  3:  Hydrostatic  Test  Configuration  [3]. 


Figure  4:  Instrumentation  of  the  Strain  Gauges  on  the  SW226  Valve  [3]. 
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3  Geometric  Configurations  of  NAB  Valves 


3.1  Introduction 

CAD  drawings  of  the  valve  geometries  were  not  available,  so  in  this  study,  it  was  required  to  first 
take  measurements  of  the  valve  dimensions  in  order  to  produce  CAD  drawings  of  the  valves. 
Due  to  the  complexity  of  the  valve  geometries,  it  was  decided  that  measurements  of  the  valve 
dimensions  be  carried  out  using  a  Coordinate  Measuring  Machine  (CMM),  rather  than  hand  tools. 
The  CMM  machine  at  the  Advanced  Manufacturing  Group  (AMG)  of  Dalhousie  University  was 
selected  for  this  purpose,  due  to  the  capabilities,  proximity  and  availability  of  the  facility  to  the 
project  team.  The  machine  is  a  Mitutoyo  BHN305  model,  which  is  a  3  axis  computer  numerical 
controlled  (CNC)  Coordinate  Measuring  Machine  (CMM)  that  is  used  to  measure  part  features 
such  as  holes,  cones,  arcs,  spheres,  surfaces,  etc.  Features  that  cannot  be  measured 
with  conventional  hand  tools  are  easily  measured  to  a  high  degree  of  precision  with  the  CMM. 
Figure  5  shows  a  picture  of  the  CMM.  It  has  the  following  capabilities 
(http://nscadcam.dal.ca/cmm  info.htm): 


Work  piece  Capacity: 

Measuring  range  (1  x  w  x  h) 

315  X  490  X  290  mm 

Working  area 

529  X  920  mm 

Maximum  work  piece  height 

365  mm 

Maximum  work  piece  mass 

150  kg 

Figure  6  shows  pictures  of  how  measurements  typical  valve  components  are  made  with  the 
CMM. 

The  CMM  directly  interfaces  with  CAD  packages,  such  as  Pro  Engineer,  to  produce  3-D 
geometries  of  the  objects.  The  CAD  software  package  will  then  provide  IGES  or  STEP  files  that 
can  be  easily  accessed  by  finite  element  packages.  CAD  drawings  of  the  valves  were  developed 
in  Pro  Engineer,  based  on  measurements  taken  with  the  CMM  machine.  Details  are  provided  in 
the  following  sections. 

3.2  SW403  Valve 

Figure  7(a)  shows  e-drawing  of  the  SW403  valve  assembly,  comprising  of  the  valve  body  (grey), 
ball  (green),  bonnet  or  valve  cover  (red),  the  spacer  (blue)  and  the  bolts  (yellow).  Only  half  of  the 
valve  assembly  is  shown.  However,  it  should  be  pointed  out  that  the  valve  was  not  completely 
symmetrical  about  the  axial  plane  that  passes  through  the  centre  of  the  valve,  as  can  be  seen  from 
the  section  view  of  the  valve  in  Figure  8(a).  Thus,  a  full  model  was  used  in  this  study  as  opposed 
to  a  symmetric  model  used  in  the  previous  study  [1],  as  discussed  later  in  Chapter  4.  In  addition, 
the  wall  thickness  of  the  thinner  portion  of  the  valve  body  was  reduced  once  by  material  removal 
to  simulate  corrosion.  Figure  7(b)  shows  e-drawing  of  the  SW403  valve  body  with  the  thickness 
reduction  region  and  Figure  8(b)  shows  the  section  view  of  the  valve  after  the  thickness 
reduction.  Detailed  dimensions  and  views  of  the  valve  assemble  are  provided  in  Annex  A. 
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3.3  SW226  Valve 


Figure  9(a)  shows  e-drawings  of  the  SW226  valve  assembly,  comprising  of  the  valve  body  (grey), 
ball  (green),  bonnet  or  valve  cover  (cyan)  the  spacer  (blue)  and  the  bolts  (yellow).  Only  half  of 
the  valve  assembly  is  shown.  However,  it  should  be  pointed  out  again  that  the  valve  was  not 
completely  symmetrical  about  the  axial  plane  that  passes  through  the  centre  of  the  valve,  as  can 
be  seen  from  the  section  view  of  the  valve  in  Figure  10(a).  Thus,  a  full  model  was  again  used  in 
this  study  as  opposed  to  a  symmetric  model  used  in  the  previous  study  [1],  as  discussed  later  in 
Chapter  4.  In  addition,  the  wall  thickness  of  the  thinner  portion  of  the  valve  body  was  gradually 
reduced  three  times  by  material  removal  to  simulate  corrosion.  Figure  9(b)  shows  e-drawing  of 
the  SW226  valve  body  with  the  thickness  reduction  region  and  Figure  10(b)  shows  the  section 
view  of  the  valve  after  the  first  thickness  reduction.  It  should  be  noted  that  the  overall 
dimensions  of  the  SW226  valve  were  larger  than  those  of  the  SW403,  but  both  valves  have 
similar  wall  thickness  around  the  valve  mid-body.  Detailed  dimensions  and  views  of  the  valve 
assemble  after  thickness  reductions  are  provided  in  Annex  A.  Details  of  the  initial  valves  can  be 
found  in  the  Appendix  A  of  [1]. 


DRDC  Atlantic  CR  2008-105 


Figure  5:  Picture  of  CMM  at  the  AMG,  Dalhousie  University. 
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Figure  6:  Picture  Showing  How  Valve  Dimensions  are  Measurements  on  the  CMM. 
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(b) 


Figure  7:  (a)  SW403  Valve  Assembly,  (b)  SW403  Valve  Body  Showing  the  Thickness  Reduction 

Region. 
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Figure  8:  (a)  Section  View  of  the  SW403  Valve  Showing  Uneven  Thickness  of  Valve  Body,  (b) 
Section  View  of  the  SW403  Valve  after  the  First  Thickness  Reduction. 
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(b) 


Figure  9:  (a)  SW226  Valve  Assembly,  (b)  SW226  Valve  Body  Showing  the  Thickness  Reduction 

Region. 
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(a) 


Figure  10:  (a)  Section  View  of  the  SW226  Valve  Showing  Uneven  Thickness  of  Valve  Body,  (b) 
Section  View  of  the  SW226  Valve  after  the  First  Thickness  Reduction. 
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4  Finite  Element  Models  of  NAB  Valves 


4.1  Finite  Element  Approach 

The  finite  element  modeling  and  analysis  approach  adopted  in  the  study  was  governed  by  the 
need  to  adequately  model  the  complex  geometric  configurations  of  the  valves  assemblies,  as  well 
as  the  nonlinear  static  behaviours  of  the  systems  under  the  intense  hydrostatic  loads.  The 
approach  consisted  of  the  use  of  the  HyperMesh  [4],  LS-DYNA  [5]  finite  element  software  suites. 
HyperMesh  was  used  to  generate  the  finite  element  models.  This  software  has  capabilities  for 
directly  importing  solid  models  from  CAD  software  such  as  Pro  Engineer  [6]  that  was  used  to 
model  the  geometric  configurations  of  the  system,  and  has  very  robust  solid  finite  element 
modeling  features.  On  the  other  hand,  the  LS-DYNA  suite  was  used  to  perform  the  static  elastic- 
plastic  responses  of  the  valve  systems.  The  software  has  a  wide  range  of  elements,  explicit  and 
implicit  solvers  for  performing  static  and  dynamic  analysis,  and  capabilities  for  handling  both 
material  and  geometric  nonlinearities.  Details  of  the  finite  element  models  and  analysis  results 
are  provided  in  the  following  sections. 

4.2  Finite  Element  Meshes 

As  discussed  in  Chapter  3,  the  geometric  configurations  of  the  valves  assemblies  were  defined 
and  modeled  in  the  Pro  Engineer  CAD  software  package.  The  geometries  were  then  exported  as 
STEP  assembly  files  for  import  into  HyperMesh. 

Due  to  the  complexity  of  the  valve  geometry,  and  in  order  to  accurately  model  its  geometry,  it 
was  decided  to  use  tetrahedral  finite  elements  to  model  the  valve  configuration.  Besides, 
Hypermesh  provides  a  convenient  auto  mesh  capability  to  generate  solid  finite  element  meshes 
for  complex  geometries,  which  is  only  available  for  meshing  with  tetrahedral  elements. 

The  tetrahedral  elements  do  not  provide  the  same  level  of  accuracy  as  the  hexahedral  (brick) 
elements.  It  is  usually  necessary  to  use  finer  tetrahedral  element  meshes  in  order  to  obtain 
reasonably  accurate  results.  To  this  end  investigations  were  performed  to  determine  the  level  of 
tetrahedral  finite  element  mesh  refinement  required  to  produce  the  same  level  of  accuracy  as  the 
hexahedral  elements.  These  preliminary  investigations  enabled  the  determination  of  the  level  of 
mesh  refinement  required  to  provide  reasonable  balance  between  accuracy  and  computational 
resources  (storage  space  and  run  times).  The  meshes  used  in  this  study  were  found  to  provide 
smooth  contours  of  displacements  and  stresses/strains. 

4.3  Mesh  Refinement 

In  the  present  study,  the  previous  finite  element  models  of  the  valves  were  refined,  in  order  to 
remove  some  of  the  simplifications  made  in  the  previous  study  [1].  As  mentioned  in  Chapter  3, 
in  the  previous  study,  symmetry  about  an  axial  plane  that  passes  through  the  centre  of  the  valve 
was  invoked  to  reduce  the  sizes  of  the  models,  in  anticipation  of  the  large  transient  nonlinear 
elastic-plastic  analyses.  The  thinner  portion  of  the  valve  was  modeled  as  this  would  govern  the 
valve  behaviour.  However,  in  the  present  study,  the  wall  thickness  of  the  thinner  portion  of  the 
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valve  body  was  gradually  reduced,  making  the  valve  increasingly  unsymmetrical  about  the  axial 
plane.  Thus,  the  modeling  of  the  full  valve  was  deemed  to  be  more  appropriate  for  the  present 
study. 

In  addition,  in  the  previous  study,  the  bonnet  was  assumed  to  be  connected  to  the  valve  body  at 
the  interface.  This  reasonable  assumption  was  made  in  order  to  simplify  the  model,  but  it  has  the 
effect  of  making  the  bonnet-valve  interface  somewhat  stiffer  than  the  physical  problem.  A  more 
accurate  representation  would  be  to  model  the  holes  and  the  bolts  in  detail.  In  the  present  study, 
the  model  was  refined  to  include  detailed  geometry  of  the  bonnet  and  the  bolt  assembly  to 
accurately  model  the  strain  in  the  bonnet.  Figure  11  shows  the  section  view  of  the  bonnet-bolt 
interface  of  a  typical  model.  An  appropriate  contact  penalty  stiffness  scale  factor  on  avoiding 
nodal  penetration  was  assigned  to  the  bonnet-bolt  interface  for  each  model. 

Figure  12  shows  the  complete  view  of  the  SW403  valve  model  with  initial  wall  thickness.  The 
wall  thickness  of  the  thinner  portion  of  this  valve  was  reduced  once,  resulting  in  a  total  of  two 
cases  of  experimental  testing  and  analyses.  Figure  13  shows  the  complete  view  of  the  SW226 
valve  model  with  initial  wall  thickness.  The  wall  thickness  of  the  thinner  portion  of  this  valve 
was  reduced  three  times,  resulting  in  a  total  of  four  cases  of  experimental  testing  and  analyses. 
The  model  statistics  are  provided  in  Table  1. 

4.4  Material  Properties 

The  valve  are  made  of  nickel  aluminum  bronze  (NAB)  material  composed  approximately  of  9- 
10%  aluminum,  5%  nickel,  5%  iron,  and  the  balance  being  copper  [2],  while  the  bolts  and  the 
valve  end  cover  plates  are  made  of  44W  steel.  The  valve  materials  can  be  assumed  to  be  both 
homogeneous  and  isotropic,  which  means  they  have  the  same  physical  and  mechanical  properties 
throughout  their  volumes  and  have  these  same  properties  in  all  directions,  such  that  the  stress- 
strain  relationship  obtained  using  standard  test  in  uni-axis  specimen  can  be  used  in  3-d  numerical 
simulation.  The  mechanical  properties  of  these  materials,  obtained  from  Refs  [7]  and  [8],  are 
summarized  in  Table  2.  All  materials  are  assumed  to  behave  in  an  elastic-plastic  manner,  and  are 
modeled  with  Material  Type  24  in  LS-DYNA  [5],  which  treats  the  stress-strain  relationship  by  a 
bilinear  stress-strain  curve  by  defining  the  yield  stress  and  the  tangent  (plastic)  modulus. 

4.5  Boundary  Conditions 

Figure  14  shows  the  boundary  conditions  used  for  a  typical  model.  The  pretension  in  the  bolts 
was  assumed  to  keep  the  top  of  the  bolts  from  moving.  Thus,  the  boundary  conditions  applied 
were: 


a.  u  =  v  =  w  =  9^=  9y  =  9z  =  0  top  of  the  bolts  and  the  bolt  contact  faces  on  the 
valve  end  cover  plates 

where  u,  v,  w  are  the  translational  displacements  in  the  global  x,  y,  and  z-directions,  respectively, 
and  9x,  9y,  9^  are  the  rotations  about  the  respective  global  axial  directions,  respectively. 
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4.6  Hydrostatic  Loading 


The  hydrostatic  load  was  applied  as  a  uniform  pressure  on  all  surfaces  inside  the  valve.  The  load 
was  incrementally  applied  from  zero  to  just  over  5000  psi  (34.47  MPa),  as  was  done  in  the  DRDC 
hydrostatic  tests. 


Table  1:  Statistics  of  the  SW403  and  SW226  Valve  Finite  Element  Models. 


Model 

Thin  Wall 
Thickness 
(mm) 

Number 
of  Nodes 

Number  of 
Tetrahedral 
Elements 

SW403 

Initial 

Thickness 

10.2 

31,243 

138,515 

Thickness 

Reduction 

5.1 

31,006 

137,180 

SW226 

Initial 

Thickness 

10.4 

56,070 

254,466 

Thickness 

Reduction 

5.5 

55,831 

253,103 

2nd 

Thickness 

Reduction 

3.5 

55,970 

253,687 

3'^  Thickness 
Reduction 

2.0 

56,186 

254,813 

Table  2:  Material  Properties. 


Property/Constant 

NAB  [7] 

44W  Steel  [8] 

Young’s  Modulus,  msi  (GPa) 

18.13  (125) 

29.7  (205) 

Yield  Stress,  ksi  (MPa) 

36.26  (250) 

43.9  (303) 

Ultimate  Strength,  ksi  (MPa) 

85  (586) 

Plastic  Modulus,  msi  (GPa) 

0.327  (2.26) 

2.97x10'^  (0.002) 

Poisson’s  Ratio 

0.32 

0.3 

Density,  Ibs^/in"^  (kg/m3) 

7.147x10'^  (7640) 

7.296x10'^  (7800) 

Elongation  (%) 

15-20 
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Figure  11:  Section  View  of  the  Bonnet-Bolt  Interface  of  a  Typical  Model. 


Figure  12:  Finite  Element  Mesh  used  for  the  SW403  Valve  with  Initial  Wall  Thickness. 
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Figure  13:  Finite  Element  Mesh  used  for  the  SW226  Valve  with  Initial  Wall  Thickness. 


Figure  14:  Boundary  Conditions  for  a  Typical  Model. 
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5  Response  of  NAB  Valves  to  Hydrostatic  Load 


The  responses  of  the  SW403  and  SW226  valves  to  the  hydrostatic  load  were  analyzed  using  the 
LS-DYNA  solver.  Linear  elastic  analysis  showed  that  the  maximum  stress  exceeded  the  NAB 
yield  stress.  Hence  a  nonlinear  elastic  plastic  analysis  was  performed  using  the  LS-DYNA 
implicit  solver.  The  load  was  applied  incrementally  in  automatically  adjusted  steps.  The  results 
of  the  analyses  are  discussed  below. 

The  resultant  displacement,  von  Mises  stress  and  effective  strain  fields  from  the  finite  element 
analysis  of  each  model  are  plotted  in  the  following  sections.  The  resultant  displacement  (R)  is  the 
magnitude  of  the  vector  sum  of  the  displacements  in  the  global  x,  y,  and  z-directions: 

R  =  ^Ju^+v^+w\  (1) 


where  u,  v  and  w  are  the  magnitudes  of  the  translational  displacement  in  the  global  x,  y,  and  z- 
directions. 

The  von  Mises  stress  ( )  is  a  scalar  stress  value  that  can  be  computed  from  the  stress  tensor. 

It  is  used  as  a  criterion  in  determining  the  onset  of  yielding  in  ductile  materials  such  as  metals. 
The  criterion  states  that  a  material  starts  yielding  when  its  von  Mises  stress  reaches  the  yield 
stress  of  the  material. 


^VM 


1 

7T 


(2) 


where  ,  (T2  and  are  the  principal  stresses. 

Similarly,  the  effective  strain  ( )  is  the  same  criterion  computed  from  the  strain  tensor: 


^VM 


(3) 


where  6*^ ,  6*2  and  6^3  are  the  principal  strains. 


5.1  SW403  Valve 

The  wall  thickness  of  the  thinner  portion  of  this  valve  was  reduced  once,  resulting  in  a  total  of 
two  cases  of  experimental  testing  and  analyses. 
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5.1.1  Initial  Wall  Thickness  of  10.2  mm 


As  indicated  in  Table  1,  the  initial  wall  thickness  of  the  thinner  portion  of  the  SW403  valve  was 

10.2  mm.  Figure  15  to  Figure  17  compare  the  contours  of  the  displacement,  von  Mises  stress  and 
effective  strain  fields,  respectively,  of  the  full  model  of  the  valve  in  the  present  study  to  the  half 
model  in  the  previous  study  [1],  at  8000  psi  (55.16  MPa)  hydrostatic  load.  Figure  18  shows  the 
von  Mises  stress  contours  of  the  yielded  areas  only.  It  can  be  seen  from  these  figures  that  the 
maximum  displacements,  stresses  and  strains  generally  occurred  in  the  centre  of  the  valve  body. 
It  can  also  be  observed  that  the  maximum  displacement,  stress  and  strain  in  the  full  model  were 
similar  to  those  in  the  half  model.  However,  the  differences  in  the  results  can  be  expected  to 
increase  as  the  wall  thickness  is  reduced  and  the  valve  becomes  more  unsymmetrical.  In  addition, 
the  maximum  load  under  which  the  strains  of  the  NAB  material  remain  in  the  elastic  region  is 
approximately  5500  psi,  as  determined  from  the  finite  element  analysis. 

Figure  19  to  Figure  21  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  the  strain  gauge  locations  to  the  half  model  and  the  experimental  results.  It  can 
be  seen  from  these  figures  that  the  load-strain  histories  from  the  models  were  comparable  to  the 
experimental  results,  except  for  the  slight  deviation  of  the  linear  behaviour  observed  in  the  neck 
strain  in  the  experiment.  It  can  also  be  observed  that  the  responses  from  the  half  model  were  in 
general  stiffer  than  the  full  model.  This  observation  is  possibly  due  to  the  fact  that  the  modeling 
of  the  bonnet- valve  interface  in  the  previous  study  was  more  constrained  than  that  in  the  present 
study. 

Figure  22  to  Figure  24  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 
be  seen  from  these  figures  that  all  three  strain  responses  remained  linear  elastic  for  the  load  range 
considered  (up  to  6000  psi).  In  addition,  the  finite  element  results  in  the  body  were  reasonably 
good  as  the  experimental  results  were  within  the  range  of  finite  element  results  near  the  strain 
gauges.  However,  the  circumferential  strain  in  the  neck  was  underestimated. 

5.1.2  After  the  First  Wall  Thickness  Reduction  to  5.1  mm 

As  indicated  in  Table  1,  the  wall  thickness  of  the  thinner  portion  of  the  SW403  valve  was  5.1  mm 
after  the  first  reduction.  Figure  25  to  Figure  27  show  the  contours  of  the  displacement,  von  Mises 
stress  and  effective  strain  fields,  respectively,  of  the  model  of  the  valve,  at  5306  psi  (36.58  MPa) 
hydrostatic  load.  Figure  28  shows  the  von  Mises  stress  contours  of  the  yielded  areas  only.  It  can 
be  seen  from  these  figures  that  the  maximum  displacements,  stresses  and  strains  again  generally 
occurred  in  the  centre  of  the  valve  body.  It  can  also  be  observed  that  the  maximum  displacement, 
stress  and  strain  in  this  model  were  greater  than  those  in  the  previous  SW403  valve  with  initial 
wall  thickness  at  similar  load,  as  expected  due  to  the  reduced  wall  thickness.  In  addition,  the 
maximum  load  under  which  the  strains  of  the  NAB  material  remain  in  the  elastic  region  decreases 
to  approximately  3200  psi,  as  determined  from  the  finite  element  analysis. 

Figure  29  to  Figure  31  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 


DRDC  Atlantic  CR  2008-105 


21 


be  seen  from  these  figures  that  the  experimental  strain  responses  showed  signs  of  yielding  in  the 
body  for  the  load  range  considered  (up  to  6000  psi).  In  addition,  within  the  linear  elastic  region, 
the  wall  thickness  reduction  caused  an  increase  of  about  34%  in  the  axial  strain  in  the  body,  an 
increase  of  about  98%  in  the  circumferential  strain  in  the  body  and  an  increase  of  about  36%  in 
the  circumferential  strain  in  the  neck  over  the  strains  in  the  valve  with  initial  wall  thickness,  as 
shown  in  Table  3,  which  summarizes  the  average  experimental  strains  at  1500  psi  hydrostatic 
load  of  each  valve.  Furthermore,  the  finite  element  results  for  the  axial  strain  in  the  body  were 
reasonably  good  as  the  experimental  results  were  within  the  range  of  finite  element  results  near 
the  strain  gauges.  However,  the  circumferential  strain  in  both  the  body  and  the  neck  were 
underestimated,  especially  in  the  body.  It  is  expected  that  further  model  refinement,  such  as  the 
use  of  higher  order  tetrahedral  elements,  may  be  required  to  improve  the  circumferential  strain 
responses.  This  could  not  be  carried  out  within  the  time  frame  of  this  project.  However,  the 
effect  of  such  model  refinement  on  the  circumferential  strains  should  be  investigated  in  the  future. 

5.2  SW226  Valve 

The  wall  thickness  of  the  thinner  portion  of  this  valve  was  reduced  three  times,  resulting  in  a  total 
of  four  cases  of  experimental  testing  and  analyses. 

5.2.1  Initial  Wall  Thickness  of  10.4  mm 

As  indicated  in  Table  1,  the  initial  wall  thickness  of  the  thinner  portion  of  the  SW226  valve  was 
10.4  mm.  Figure  32  to  Figure  34  show  the  contours  of  the  displacement,  von  Mises  stress  and 
effective  strain  fields,  respectively,  of  the  model  of  the  valve,  at  5323  psi  (36.70  MPa)  hydrostatic 
load.  Figure  35  shows  the  von  Mises  stress  contours  of  the  yielded  areas  only.  It  can  be  seen 
from  these  figures  that  the  maximum  displacements,  stresses  and  strains  again  generally  occurred 
in  the  centre  of  the  valve  body.  It  can  also  be  observed  that  the  maximum  displacement,  stress 
and  strain  in  this  valve  were  greater  than  those  in  the  smaller  SW403  valve  at  similar  load.  An 
immediate  implication  of  this  observation  is  that  corrosion  of  the  valve  material  around  the  mid¬ 
body  in  SW226  could  have  more  pronounced  effects  on  the  structural  integrity  than  in  SW403.  In 
addition,  the  maximum  load  under  which  the  strains  of  the  NAB  material  remain  in  the  elastic 
region  is  approximately  3600  psi,  as  determined  from  the  finite  element  analysis. 

Figure  36  to  Figure  38  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 
be  seen  from  these  figures  that  all  three  experimental  strain  responses  remained  linear  elastic  for 
the  load  range  considered  (up  to  6000  psi).  Comparing  the  strains  of  the  two  valves  with  initial 
wall  thickness  in  Table  3,  it  can  again  be  seen  that  the  smaller  SW403  valve  was  stiffer  and 
responded  with  lower  strains  than  the  larger  SW226  valve.  In  addition,  the  finite  element  results 
for  the  axial  strain  in  the  body  were  reasonably  good  as  the  experimental  results  were  within  the 
range  of  finite  element  results  near  the  strain  gauges.  However,  the  circumferential  strains  were 
again  underestimated,  suggesting  further  model  refinement  may  be  required. 
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5.2.2  After  the  First  Wall  Thickness  Reduction  to  5.5  mm 


As  indicated  in  Table  1,  the  wall  thickness  of  the  thinner  portion  of  the  SW226  valve  was  5.5  mm 
after  the  first  reduction.  Figure  39  to  Figure  41  show  the  contours  of  the  displacement,  von  Mises 
stress  and  effective  strain  fields,  respectively,  of  the  model  of  the  valve,  at  5278  psi  (36.39  MPa) 
hydrostatic  load.  Figure  42  shows  the  von  Mises  stress  contours  of  the  yielded  areas  only.  It  can 
be  seen  from  these  figures  that  the  maximum  displacements,  stresses  and  strains  again  generally 
occurred  in  the  centre  of  the  valve  body.  It  can  also  be  observed  that  the  maximum  displacement, 
stress  and  strain  in  this  model  were  greater  than  those  in  the  previous  SW226  valve  with  initial 
wall  thickness  at  similar  load,  as  expected  due  to  the  reduced  wall  thickness.  In  addition,  the 
maximum  load  under  which  the  strains  of  the  NAB  material  remain  in  the  elastic  region  decreases 
to  approximately  3300  psi,  as  determined  from  the  finite  element  analysis. 

Figure  43  to  Figure  45  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 
be  seen  from  these  figures  that  the  both  the  experimental  and  finite  element  strain  responses 
showed  signs  of  yielding  in  the  body  for  the  load  range  considered  (up  to  6000  psi).  In  addition, 
within  the  linear  elastic  region,  the  wall  thickness  reduction  caused  an  increase  of  about  44%  in 
the  axial  strain  in  the  body,  an  increase  of  about  37%  in  the  circumferential  strain  in  the  body  and 
a  decrease  of  about  1 1%  in  the  circumferential  strain  in  the  neck  over  the  strains  in  the  valve  with 
initial  wall  thickness,  as  shown  in  Table  3.  Comparing  the  strains  of  the  two  valves  after  the  first 
wall  thickness  reduction  in  Table  3,  it  is  interesting  to  note  that  while  the  experimental  axial  strain 
in  the  body  and  the  circumferential  strain  in  the  neck  of  the  smaller  and  stiffer  SW403  valve  were 
lower  than  this  valve,  the  circumferential  strain  in  the  body  was  higher,  which  had  similar  wall 
thickness  (5.1  mm  for  SW403  and  5.5  mm  for  SW226).  This  indicates  that  corrosion  damage  has 
a  more  significant  effect  on  the  circumferential  strain  of  the  SW403  valve  body  than  of  the 
SW226  valve,  causing  a  higher  response  from  the  SW403  valve  when  a  stiffer  response  was 
expected.  Furthermore,  the  finite  element  results  underestimated  all  three  experimental  strain 
responses,  especially  for  the  circumferential  strain  in  the  body.  Nonetheless,  the  finite  element 
results  appeared  to  be  able  to  predict  the  onset  of  yielding  in  the  circumferential  strain  in  the 
body. 

5.2.3  After  the  Second  Wall  Thickness  Reduction  to  3.5  mm 

As  indicated  in  Table  1,  the  wall  thickness  of  the  thinner  portion  of  the  SW226  valve  was  3.5  mm 
after  the  second  reduction.  Figure  46  to  Figure  48  show  the  contours  of  the  displacement,  von 
Mises  stress  and  effective  strain  fields,  respectively,  of  the  model  of  the  valve,  at  5323  psi  (36.70 
MPa)  hydrostatic  load.  Figure  49  shows  the  von  Mises  stress  contours  of  the  yielded  areas  only. 
It  can  be  seen  from  these  figures  that  the  maximum  displacements,  stresses  and  strains  again 
generally  occurred  in  the  centre  of  the  valve  body.  It  can  also  be  observed  that  the  maximum 
displacement,  stress  and  strain  in  this  model  were  greater  than  those  in  the  previous  SW226  cases 
at  similar  load,  as  expected  due  to  the  reduced  wall  thickness.  In  addition,  the  maximum  load 
under  which  the  strains  of  the  NAB  material  remain  in  the  elastic  region  further  decreases  to 
approximately  2200  psi,  as  determined  from  the  finite  element  analysis. 
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Figure  50  to  Figure  52  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 
be  seen  from  these  figures  that  the  both  the  experimental  and  finite  element  strain  responses  again 
showed  signs  of  yielding  in  the  body  for  the  load  range  considered  (up  to  4000  psi),  which  was 
lower  than  that  of  the  previous  SW226  cases.  In  addition,  within  the  linear  elastic  region,  the 
wall  thickness  reduction  caused  an  increase  of  about  91%  in  the  axial  strain  in  the  body,  an 
increase  of  about  223%  in  the  circumferential  strain  in  the  body  and  a  decrease  of  about  19%  in 
the  circumferential  strain  in  the  neck  over  the  strains  in  the  valve  with  initial  wall  thickness,  as 
shown  in  Table  3.  Thus,  it  can  be  seen  that  further  corrosion  damage  had  the  most  significant 
effect  on  the  circumferential  strain  of  the  valve  body,  causing  it  to  dominate  the  valve  response  to 
the  hydrostatic  loading.  In  addition,  the  corrosion  damage  effect  on  the  neck  strain  remained  less 
significant  compared  to  the  body  strains.  Furthermore,  the  finite  element  results  for  the  axial 
strain  in  the  body  and  the  circumferential  strain  in  the  neck  were  reasonably  good  as  the 
experimental  results  were  within  the  range  of  finite  element  results  near  the  strain  gauges. 
However,  the  circumferential  strain  in  the  body  was  again  underestimated.  Nonetheless,  the  finite 
element  results  appeared  to  be  able  to  predict  the  onset  of  yielding  in  the  circumferential  strain  in 
the  body. 

5.2.4  After  the  Third  Wall  Thickness  Reduction  to  2.0  mm 

As  indicated  in  Table  1,  the  wall  thickness  of  the  thinner  portion  of  the  SW226  valve  was  2.0  mm 
after  the  third  reduction.  Figure  53  to  Figure  55  show  the  contours  of  the  displacement,  von 
Mises  stress  and  effective  strain  fields,  respectively,  of  the  model  of  the  valve,  at  5337  psi  (36.80 
MPa)  hydrostatic  load.  Figure  56  shows  the  von  Mises  stress  contours  of  the  yielded  areas  only. 
It  can  be  seen  from  these  figures  that  the  maximum  displacements,  stresses  and  strains  again 
generally  occurred  in  the  centre  of  the  valve  body.  It  can  also  be  observed  that  the  maximum 
displacement,  stress  and  strain  in  this  model  were  greater  than  those  in  the  previous  SW226  cases 
at  similar  load,  as  expected  due  to  the  reduced  wall  thickness.  In  addition,  the  maximum  load 
under  which  the  strains  of  the  NAB  material  remain  in  the  elastic  region  further  decreases  to 
approximately  1400  psi,  as  determined  from  the  finite  element  analysis. 

Figure  57  to  Figure  59  compare  the  load-axial  strain  history  in  the  body,  the  load-circumferential 
strain  history  in  the  body  and  the  load-circumferential  strain  history  in  the  neck,  respectively,  of 
the  full  model  at  selected  locations  near  the  strain  gauges  (SG)  to  the  experimental  results.  It  can 
be  seen  from  these  figures  that  all  three  experimental  strain  responses  remained  linear  elastic  for 
the  load  range  considered  (up  to  2500  psi),  which  was  lower  than  that  of  the  previous  SW226 
cases.  In  addition,  within  the  linear  elastic  region,  the  wall  thickness  reduction  caused  an  increase 
of  about  140%  in  the  axial  strain  in  the  body,  an  increase  of  about  427%  in  the  circumferential 
strain  in  the  body  and  a  decrease  of  about  25%  in  the  circumferential  strain  in  the  neck  over  the 
strains  in  the  valve  with  initial  wall  thickness,  as  shown  in  Table  1.  Thus,  it  can  again  be  seen 
that  further  corrosion  damage  had  the  most  significant  effect  on  the  circumferential  strain  of  the 
valve  body,  causing  it  to  dominate  the  valve  response  to  the  hydrostatic  loading.  In  addition,  the 
corrosion  damage  effect  on  the  neck  strain  again  remained  less  significant  compared  to  the  body 
strains.  Furthermore,  the  finite  element  results  for  the  axial  strain  in  the  body  and  the 
circumferential  strain  in  the  neck  were  reasonably  good  as  the  experimental  results  were  within 
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the  range  of  finite  element  results  near  the  strain  gauges.  However,  the  circumferential  strain  in 
the  body  was  again  underestimated. 


5.3  Effect  of  Corrosion  Degradation  on  Vaive  Performance 

The  goal  of  the  present  study  is  to  evaluate  the  effect  of  corrosion  degradation  on  the  performance 
of  the  NAB  valves.  In  the  present  study,  this  can  be  expressed  in  terms  of  performance  charts, 
which  plot  the  various  strain  responses  against  the  valve  body  wall  thickness.  Unfortunately, 
there  are  not  enough  data  in  the  present  study  to  generate  performance  charts  for  the  SW403 
valve.  Nonetheless,  the  performance  charts  for  the  axial  and  the  circumferential  strain  in  the 
SW226  valve  body  at  1500  psi  (10.34  MPa)  hydrostatic  load  are  shown  in  Figure  60  and  Figure 
61,  respectively.  The  performance  of  the  strains  in  the  neck  was  deemed  insignificant  when 
compared  to  the  strains  in  the  body.  It  can  be  seen  from  the  figures  that  the  agreement  between 
the  experimental  and  finite  element  axial  strain  results  was  reasonably  good,  while  the  finite 
element  results  underestimated  the  circumferential  strains.  Thus,  further  model  refinement  would 
be  required  in  order  to  generate  the  complete  performance  charts  for  the  valves. 

The  approximate  maximum  loads  under  which  the  strains  of  the  NAB  material  in  the  valves 
remain  in  the  elastic  region,  as  determined  from  the  finite  element  analyses,  are  summarized  in 
Table  4.  It  can  be  seen  that  the  load  decreased  with  corrosion  damage,  as  expected.  The  results 
for  the  SW226  valve  are  plotted  in  Figure  62. 
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Table  3:  Comparison  of  Experimental  Strains  of  the  SW403  and  SW226  Valves  at  1500  psi  (10.34 

MPa)  Hydrostatic  Load. 


Valve 

Axial 

Body 

%  Change 
from 
Initial 

Strain 

Circum¬ 

ferential 

Body 

microstrain) 
%  Change 
from 
Initial 

Circum¬ 

ferential 

Neck 

%  Change 
from 
Initial 

SW403 

Initial 
Thickness 
(10.2  mm) 

237 

311 

112 

Thickness 
Reduction 
(5.1  mm) 

317 

33.8 

617 

98.2 

152 

36.0 

SW226 

Initial 
Thickness 
(10.4  mm) 

484 

371 

217 

1^^  Thickness 
Reduction 
(5.5  mm) 

699 

44.4 

506 

36.5 

192 

-11.4 

2^^  Thickness 
Reduction 
(3.5  mm) 

925 

91.1 

1,196 

222.8 

175 

-19.3 

3^^  Thickness 
Reduction 
(2.0  mm) 

1,162 

140.0 

1,952 

426.7 

164 

-24.6 

Table  4:  Approximate  Maximum  Load  under  Which  the  Strains  of  the  NAB  Material  Remain  in 
the  Elastic  Region  as  Determined  from  the  Finite  Element  Analyses. 


Valve 

Maximum  Load  for 
Elastic  Strains  (psi) 

SW403 

Initial  Thickness  (10.2  mm) 

5,500 

1^^  Thickness  Reduction  (5.1  mm) 

3,200 

SW226 

Initial  Thickness  (10.4  mm) 

3,600 

1^^  Thickness  Reduction  (5.5  mm) 

3,300 

2"‘*  Thickness  Reduction  (3.5  mm) 

2,200 

3'^‘*  Thickness  Reduction  (2.0  mm) 

1,400 
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Fringe  Levels 


Time  =  1 

Contours  of  Resultant  Displacement  (in) 
min=3.00479e-05,  at  node#  30328 


8.495e-03 

7.648e-03 

6.802e-03 

5.955e-03 

5.109e-03  _ 

4.262e4)3  _ 

3.416e-03 

2.569e-03  _ 

1.723e4)3 

8.765e-04 

3.005e-05 


1 

1 


(a) 


Fringe  Levels 
8.702e-03 
7.832e4)3 
6.961  e-03 
6.091  e4)3 
5.221  e4)3 
4.351  e4)3 
3.481  e4]3 
2.611e4)3 
1.740e4)3 
8.702e4)4 
0.000e«00 


(b) 

Figure  15:  Displacement  Contours  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm)  at 
8000  psi  (55.16  MPa)  Hydrostatic  Load:  (a)  Full  Model,  (b)  Half  Model. 
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Fringe  Levels 


Time  =  1 

Contours  of  Effective  Stress  (v-m)  (psi) 
min=99.3611,  at  elemA'  96876 
max=36797.4, 


3.680e«04 

3.313e+04 

2.946e^04 

2.579e+04 

2.212e+04 

1.845e*04 

1.478e+04 

1.111e+04 

7.439e+B3 

3.769e+B3 

9.936e*01 


) 

] 


(a) 


"T 


(b) 

Figure  16:  von  Mises  Stress  Contours  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm) 
at  8000  psi  (55.16  MPa)  Hydrostatic  Load:  (a)  Full  Model,  (b)  Half  Model. 
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Time  =  1 

Contours  of  Effective  Strain  (in/in) 
m in=4 .82282 e -06,  at  eiem#  96876 


Fringe  Leveis 
3.398e4)3 
3.058e-03 
2.719e-03 
2.380e-03 
2.041  e-03 
1.701e4)3 
1.362e-03 
1.023e4)3 
6.834e-04 
3.441  e-04 
4.823e-06 


(a) 


Fringe  Levels 


3.400e-03 
3.061  e-03 
2.722e-03 
2.383e-03 
2.043e4)3 
1.704e-03 
1.365e-03 
1.026e-03 
6.870e-04 
3.479e-04 
8.724e-06 


(b) 

Figure  17:  Effective  Strain  Contours  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm)  at 
8000  psi  (55.16  MPa)  Hydrostatic  Load:  (a)  Full  Model,  (b)  Half  Model. 
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Time  =  1 

Contours  of  Effective  Stress  (v-m)  (psi) 
min=36260,  at  elem#  96876 


Fringe  Levels 
3.680e^04 
3.674e*04 
3.669e+04 
3.664e*04 
3.658e+04 
3.653e+04 
3.647e+04 
3.642e+04 
3.637e+04 
3.631  e+04 
3.626e^04 


1 


1 


(b> 

Figure  18:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW403  Valve  with  Initial 
Wall  Thickness  (10.2  mm)  at  8000  psi  (55.16  MPa)  Hydrostatic  Load:  (a)  Full  Model,  (b)  Half 

Model. 
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0  200  400  600  800  1000 

Axial  Strain  (microstrain) 

Figure  19:  Comparison  of  Experimental,  Full  Model  and  Half  Model  Load-Axial  Strain  History 
at  the  Body  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 


Load-Circumferential  Strain  (SW403  Body,  Initial  Thickness) 


0  200  400  600  800  1000  1200  1400 

Circumferential  Strain  (microstrain) 


Figure  20:  Comparison  of  the  Experimental,  Full  Model  and  Half  Model  Load-Circumferential 
Strain  History  at  the  Body  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 


DRDC  Atlantic  CR  2008-105 


31 


Load-Circumferential  Strain  (SW403  Neck,  Initial  Thickness) 


0  100  200  300  400  500 

Circumferential  Strain  (microstrain) 


Figure  21:  Comparison  of  the  Experimental,  Full  Model  and  Half  Model  Load-Circumferential 
Strain  History  at  the  Neck  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 


Load-Axial  Strain  (SW403  Body,  Initial  Thickness) 
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Figure  22:  Comparison  of  the  Experimental  and  FE  Load-Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 
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Load-Circumferential  Strain  (SW403  Body,  Initial  Thickness) 


0  200  400  600  800  1000  1200  1400 

Circumferential  Strain  (microstrain) 


Figure  23:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 


Load-Circumferential  Strain  (SW403  Neck,  Initial  Thickness) 


0  100  200  300  400  500 

Circumferential  Strain  (microstrain) 


Figure  24:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm). 
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Time  =  0.88425 


Fringe  Levels 


Contours  of  Resultant  Displacement  (in) 
min-0.000306577,  at  nodeiV  28330 
max=0.010301, 


1.030e4)2 

9.302e4)3 

8.302e4)3 

7.303e4)3 

6.303e4)3 

5.304e4)3 

4.304e4)3 

3.305e4)3 

2.305e4)3 

1.306e4)3 

3.066e4)4 


1 


Figure  25:  Displacement  Contours  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction 
(5.1  mm)  at  5306 psi  (36.58  MPa)  Hydrostatic  Load. 


Time  =  0.88425 

Contours  of  Effective  Stress  (v-m)  (psi) 
min=63.2596,  at  elem#  23046 
max=36962.1, 


Fringe  Levels 
3.696e«04 
3.327e+04 
2.958e^04 
2.589e«04 
2.220e+04 
1.851e+04 
1.482e+04 
1.113e*04 
7.443e+03 
3.753e+03 
6.326e^01 


Figure  26:  von  Mises  Stress  Contours  of  the  SW403  Valve  after  the  First  Wall  Thickness 
Reduction  (5.1  mm)  at  5306 psi  (36.58  MPa)  Hydrostatic  Load. 
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Time  =  0.88425 

Contours  of  Effective  Strain  (in/in) 

min=3.07052e4)6,  at  elem#  23046 


Fringe  Levels 
3.901  e4)3 
3.511e4)3 
3.121e4)3 
2.731  e4)3 
2.342e4)3 
1.952e4)3 
1.562e4)3 
1.172e4)3 
7.826e4)4 
3.928e4)4 
3.071  e4)6 


Figure  27:  Effective  Strain  Contours  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction 
(5.1  mm)  at  5306 psi  (36.58  MPa)  Hydrostatic  Load. 


Time  =  0.88425 

Contours  of  Effective  Stress  (v-m)  (psi) 
min=36260,  at  elem#  23046 


Fringe  Levels 
3.696e«04 
3.689e«04 
3.682e^04 
3.675e+04 
3.668e«04 
3.661  e«04 
3.654e«04 
3.647e+04 
3.640e+04 
3.633e+04 
3.626e«04 


) 


] 


Figure  28:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW403  Valve  after  the  First 
Wall  Thickness  Reduction  (5.1  mm)  at  5306 psi  (36.58  MPa)  Hydrostatic  Load. 
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Load-Axial  Strain  (SW403  Body,  1st  Reduction) 
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Figure  29:  Comparison  of  the  Experimental  and  FE  Load-Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction  (5.1  mm). 
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Figure  30:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction  (5.1 

mm). 
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Load-Circumferential  Strain  (SW403  Neck,  1st  Reduction) 
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Figure  31:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction  (5.1 

mm). 


Time  =  0.8871 

Contours  of  Resultant  Displacement  (in) 
min=0.000193301,  at  node#  29514 
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Figure  32:  Displacement  Contours  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm)  at 

5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Time  =  0.8871 

Contours  of  Effective  Stress  (v-m) 
min=67.0859,  at  elem#  287356 
max=37349.9, 


Fringe  Levels 
3.735e+04 
3.362e+04  _ 
2.989e+04  _ 
2.617e+04  _ 
2.244e+04  _ 
1.871e+04  _ 
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Figure  33:  von  Mises  Stress  Contours  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm) 

at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 


Time  =  0.8871 

Contours  of  Effective  Strain  (in/in) 
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Figure  34:  Effective  Strain  Contours  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm)  at 

5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Time  =  0.8871 

Contours  of  Effective  Stress  (v-m)  (psi) 
min-36260,  at  elemA'  287356 


Fringe  Levels 
3.735e+04 
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Figure  3S:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW226  Valve  with  Initial 
Wall  Thickness  (10.4  mm)  at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 


Load-Axial  Strain  (SW226  Body,  Initial  Thickness) 


Figure  36:  Comparison  of  the  Experimental  and  FE  Load- Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm). 
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Load-Circumferential  Strain  (SW226  Body,  Initial  Thickness) 
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Figure  37:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm). 
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Figure  38:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm). 
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Fringe  Levels 


Time  =  0.87968 

Contours  of  Resultant  Displacement  (in) 
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Figure  39:  Displacement  Contours  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction 
(5.5  mm)  at  5278 psi  (36.39  MPa)  Hydrostatic  Load. 
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Figure  40:  von  Mises  Stress  Contours  of  the  SW226  Valve  after  the  First  Wall  Thickness 
Reduction  (5.5  mm)  at  5278 psi  (36.39  MPa)  Hydrostatic  Load. 
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Figure  41:  Effective  Strain  Contours  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction 
(5.5  mm)  at  5278 psi  (36.39  MPa)  Hydrostatic  Load. 
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Figure  42:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW226  Valve  after  the  First 
Wall  Thickness  Reduction  (5.5  mm)  at  5278 psi  (36.39  MPa)  Hydrostatic  Load. 


42 


DRDC  Atlantic  CR  2008-105 


Load-Axial  Strain  (SW226  Body,  1st  Reduction) 
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Figure  43:  Comparison  of  the  Experimental  and  FE  Load-Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction  (5.5  mm). 
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Figure  44:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction  (5.5 

mm). 
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Load-Circumferential  Strain  (SW226  Neck,  1st  Reduction) 
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Figure  45:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction  (5.5 

mm). 
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Figure  46:  Displacement  Contours  of  the  SW226  Valve  after  the  Second  Wall  Thickness 
Reduction  (3.5  mm)  at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Figure  47:  von  Mises  Stress  Contours  of  the  SW226  Valve  after  the  Second  Wall  Thickness 
Reduction  (3.5  mm)  at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Figure  48:  Effective  Strain  Contours  of  the  SW226  Valve  after  the  Second  Wall  Thickness 
Reduction  (3.5  mm)  at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Figure  49:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW226  Valve  after  the 
Second  Wall  Thickness  Reduction  (3.5  mm)  at  5323  psi  (36. 70  MPa)  Hydrostatic  Load. 
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Figure  50:  Comparison  of  the  Experimental  and  FE  Load-Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW226  Valve  after  the  Second  Wall  Thickness  Reduction  (3.5  mm). 
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Load-Circumferential  Strain  {SW226  Body,  2nd  Reduction) 
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Figure  51:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW226  Valve  after  the  Second  Wall  Thickness  Reduction 

(3.5  mm). 
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Figure  52:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW226  Valve  after  the  Second  Wall  Thickness  Reduction 

(3.5  mm). 
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Figure  53:  Displacement  Contours  of  the  SW226  Valve  after  the  Third  Wall  Thickness  Reduction 
(2.0  mm)  at  5337 psi  (36.80  MPa)  Hydrostatic  Load. 
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Figure  54:  von  Mises  Stress  Contours  of  the  SW226  Valve  after  the  Third  Wall  Thickness 
Reduction  (2.0  mm)  at  5337 psi  (36.80  MPa)  Hydrostatic  Load. 
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Figure  55:  Effective  Strain  Contours  of  the  SW226  Valve  after  the  Third  Wall  Thickness 
Reduction  (2.0  mm)  at  5337 psi  (36.80  MPa)  Hydrostatic  Load. 
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Figure  56:  von  Mises  Stress  Contours  of  the  Yielded  Areas  of  the  SW226  Valve  after  the  Third 
Wall  Thickness  Reduction  (2.0  mm)  at  5337 psi  (36.80  MPa)  Hydrostatic  Load. 
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Figure  57:  Comparison  of  the  Experimental  and  FE  Load- Axial  Strain  History  at  Selected 
Locations  on  the  Body  of  the  SW226  Valve  after  the  Third  Wall  Thickness  Reduction  (2.0  mm). 
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Figure  58:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Body  of  the  SW226  Valve  after  the  Third  Wall  Thickness  Reduction  (2.0 

mm). 
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Load-Circumferential  Strain  {SW226  Neck,  3rd  Reduction) 


2500 

2000 

^  1500 

a 

■D 

cu 

5  1000 
500 

0 


Circumferential  Strain  (microstrain) 


Figure  59:  Comparison  of  the  Experimental  and  FE  Load-Circumferential  Strain  History  at 
Selected  Locations  on  the  Neck  of  the  SW226  Valve  after  the  Third  Wall  Thickness  Reduction  (2.0 

mm). 


SW226  Body  Axial  Strain  at  1500  psi  (10.34  MPa) 


Figure  60:  Effect  of  Corrosion  Degradation  on  the  Axial  Strain  of  the  SW226  Valve  Body. 
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SW226  Body  Circumferential  Strain  at  1500  psi  (10.34  MPa) 


Figure  61:  Effect  of  Corrosion  Degradation  on  the  Circumferential  Strain  of  the  SW226  Valve 

Body. 
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Figure  62:  Effect  of  Corrosion  Degradation  on  the  Maximum  Load  under  Which  the  Strains  of 
the  SW226  NAB  Material  Remain  in  the  Elastic  Region  as  Determined  from  the  Finite  Element 

Analyses. 
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6  Summary,  Conclusions  and  Recommendations 

6.1  Summary  and  Conclusions 

This  study  was  undertaken  in  support  of  on-going  mid-life  evaluations  of  the  Victoria  Class 
submarine  valve  systems  at  DRDC  Atlantic  and  Suffield.  In  order  for  the  valves  to  perform  their 
intended  functions,  it  is  important  that  the  integrity  (structural  and  water  tightness)  of  the  valves 
is  maintained  under  all  operational  loading  conditions  as  well  as  combat  and  accidental  loads, 
such  as  shock  and  under  water  explosions,  throughout  the  life  of  the  vessel.  The  objective  of  this 
study  was  to  numerically  investigate  the  behaviour  of  the  submarine  nickel  aluminum  bronze 
(NAB)  valves  with  corrosion  damage  under  hydrostatic  loading  condition,  and  to  validate  the 
numerical  models  using  experimental  results.  Two  valve  systems,  designated  as  SW403  and 
SW226,  were  investigated.  Corrosion  damage  was  simulated  by  gradually  removing  material 
from  the  inside  the  valve  body.  CAD  drawings  of  these  valves  were  developed  by  taking 
measurements  of  the  valve  dimensions  using  a  Coordinate  Measurement  Machine  (CMM).  Finite 
element  models  of  the  valve  systems  were  then  developed  to  simulate  the  hydrostatic  test 
configuration  developed  at  DRDC  Atlantic.  Finite  element  analysis,  involving  nonlinear  elastic 
plastic  static  analysis  was  performed  using  the  implicit  solver  in  LS-DYNA  software  system. 

Under  hydrostatic  loading,  the  experimental  testing  recorded  the  axial  and  the  circumferential 
strains  in  the  valve  body  and  the  circumferential  strain  in  the  neck  of  each  valve.  The 
experimental  testing  showed  that  the  smaller  (50  lb  (24  kg))  SW403  valve  generally  responded 
with  lower  strains  than  the  larger  (100  lb  (45  kg))  SW226  valve  at  the  same  load,  except  for  the 
circumferential  strain  in  the  valve  body  of  the  SW403  valve  after  corrosion  damage.  In  addition, 
the  strains  at  the  neck  were  much  lower  than  those  at  the  body,  especially  with  increased 
corrosion  damage.  For  the  SW403  valve,  the  circumferential  strain  in  the  body  was  dominant, 
and  all  three  strain  responses  increased  with  corrosion  damage.  For  the  SW226  valve,  the  axial 
strain  in  the  body  was  initially  dominant.  However,  although  both  the  axial  and  the 
circumferential  strains  in  the  body  increased  with  corrosion  damage,  the  circumferential  strain 
increased  in  a  faster  rate  became  dominant  at  some  point.  Meanwhile,  the  circumferential  strain 
in  the  neck  decreased  with  corrosion  damage.  However,  both  magnitude  and  the  change  in  the 
neck  strain  were  less  significant  when  compared  to  the  body  strains.  This  is  expected  since  there 
was  no  thickness  reduction  in  the  neck  region,  there  should  be  no  effect  of  corrosion  damage  on 
the  circumferential  strain  in  the  neck. 

In  the  present  study,  the  previous  finite  element  models  of  the  valves  were  refined,  in  order  to 
remove  some  of  the  simplifications  made  in  the  previous  study.  In  the  present  study,  the  wall 
thickness  of  the  thinner  portion  of  the  valve  body  was  gradually  reduced  to  simulate  corrosion 
damage,  making  the  valve  increasingly  unsymmetrical  about  the  axial  plane.  Thus,  the  modeling 
of  the  full  valve  was  deemed  to  be  more  appropriate  for  the  present  study.  In  addition,  the  model 
was  refined  to  include  detailed  geometry  of  the  bonnet  and  the  bolt  assembly  to  accurately  model 
the  strain  in  the  bonnet.  An  appropriate  contact  penalty  stiffness  scale  factor  on  avoiding  nodal 
penetration  was  assigned  to  the  bonnet-bolt  interface  for  each  model.  The  finite  element  analysis 
of  the  SW403  valve  with  initial  wall  thickness  showed  that  the  results  in  the  present  model  were 
similar  to  those  in  the  previous  model.  However,  the  differences  in  the  results  can  be  expected  to 
increase  as  the  wall  thickness  is  reduced  and  the  valve  becomes  more  unsymmetrical.  In  addition. 
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the  finite  element  analyses  showed  that  the  maximum  displacements,  stresses  and  strains 
generally  occurred  in  the  centre  of  the  valve  body,  with  the  smaller  SW403  valve  generally 
responding  with  lower  values  than  the  larger  SW226  valve  at  the  same  load.  Again,  the  strains  at 
the  neck  were  much  lower  than  those  at  the  body.  Furthermore,  the  finite  element  results  for  the 
axial  strain  in  the  body  were  in  general  reasonably  good  as  the  experimental  results  were  within 
the  range  of  finite  element  results  near  the  strain  gauges.  However,  the  circumferential  strain  in 
both  the  body  and  the  neck  were  in  general  underestimated,  especially  in  the  body.  This  indicates 
that  further  refinement  of  the  finite  element  mesh  is  required. 

Finally,  in  the  experiment,  repeated  testing  was  conducted  under  various  hydrostatic  loading 
conditions  and  residual  plastic  strains  were  observed  in  the  valves  with  reduced  wall  thickness. 
However,  in  the  comparison  between  the  experimental  and  numerical  modeling  results,  these 
residual  strains  were  not  taken  into  account.  Nonetheless,  the  experimental  data  showed  that 
there  was  little  change  in  the  linear  portion  of  the  load-strain  curves  between  repeated  testing  at 
the  same  wall  thickness,  which  is  expected.  Therefore,  the  comparison  between  the  experimental 
and  numerical  modeling  results  should  be  valid  for  the  linear  portion  of  the  load-strain  curves  at 
least,  even  without  accounting  for  the  repeated  testing  and  residual  strains  in  the  numerical 
modeling.  On  the  other  hand,  in  two  cases  (the  circumferential  strain  at  the  SW226  body  after  the 
first  and  the  second  wall  thickness  reduction),  the  yield  stress  appeared  to  have  increased  by  the 
repeated  testing,  which  is  also  expected  but  could  not  be  captured  by  the  present  numerical 
modeling.  One  recommendation  on  how  to  improve  the  test  method  to  obtain  the  best  data  for 
comparison  purpose  is  to  keep  the  load  within  the  linear  elastic  region  of  the  valve  material  in 
order  to  not  introduce  residual  plastic  strains  in  the  valves. 

6.2  Recommendations 

The  following  are  recommendations  for  future  work: 

1.  The  effects  of  issues  such  as  mesh  size,  element  formulation,  element  order  and  material 
properties  should  be  investigated.  For  example,  the  underestimation  of  many  of  the 
experimental  strains  by  the  FE  models  suggests  that  the  models  were  too  stiff,  such  that  the 
use  of  a  finer  mesh,  an  element  formulation  with  nodal  rotation  or  higher  order  elements 
might  improve  the  FE  results.  Addressing  these  issues  would  improve  the  accuracy  of  the  FE 
models  so  that  they  more  closely  match  experimental  measurements  and  further  increase 
confidence  in  their  reliability. 

2.  An  assessment  of  other  valve  systems,  especially  the  ones  for  which  hydrostatic  testing  and 
UNDEX  test  results  are  available,  should  be  performed.  Numerical  models  of  these  valves 
should  also  be  performed  to  further  validate  the  numerical  approach.  Once  properly 
validated,  the  numerical  models  can  be  used  to  generate  performance  charts  for  existing  and 
new  valve  designs. 

3.  Experimentally-validated  FE  models  for  NAB  valves  could  be  extended  to  evaluate  valve 
survivability  under  conditions  beyond  hydrostatic  and  SVTV  UNDEX  tests.  For  example, 
the  SVTV  experiments  were  designed  based  on  UNDEX  loading  to  hull-mounted  valves, 
however,  FE  models  could  be  used  to  evaluate  shock  loading  to  valves  mounted  on  pipes. 
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4.  More  data  points  should  be  generated  to  enhance  corrosion  degradation  performance  curves 
for  the  SW403  valve.  In  addition,  the  curves  should  be  further  refined  after  the 
improvements  in  the  circumferential  strain  accuracy. 
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Annex  A  CAD  Drawings  of  NAB  Valves  after  Wall 
Thickness  Reductions 


(a)  (b) 

Figure  A-63;  Section  View  of  the  SW403  Valve  with  Initial  Wall  Thickness  (10.2  mm):  (a) 

Vertical  View,  (b)  Horizontal  View. 


(a)  (b) 

Figure  A-64:  Section  View  of  the  SW403  Valve  after  the  First  Wall  Thickness  Reduction  (5.1 
mm):  (a)  Vertical  View,  (b)  Horizontal  View. 
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(a)  (b) 

Figure  A-65:  Section  View  of  the  SW226  Valve  with  Initial  Wall  Thickness  (10.4  mm):  (a) 

Vertical  View,  (b)  Horizontal  View. 


(a)  (b) 

Figure  A-66:  Section  View  of  the  SW226  Valve  after  the  First  Wall  Thickness  Reduction  (5.5 
mm):  (a)  Vertical  View,  (b)  Horizontal  View. 
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(a)  (b) 

Figure  A-67:  Section  View  of  the  SW226  Valve  after  the  Second  Wall  Thickness  Reduction  (3.5 

mm):  (a)  Vertical  View,  (b)  Horizontal  View. 


(a)  (b) 

Figure  A-68:  Section  View  of  the  SW226  Valve  after  the  Third  Wall  Thickness  Reduction  (2.0 
mm):  (a)  Vertical  View,  (b)  Horizontal  View. 
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